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Bare-state population inversion is demonstrated in a two-level system with all dipole matrix 
elements nonzero. A laser field is resonantly driving the sample whereas a second weaker and lower 
frequency coherent field additionally pumps it near resonance with the dynamically-Stark-splitted 
states. Due to existence of differing permanent dipole moments in the excited and ground bare 
states, quantum coherences among the involved dressed-states are induced leading to inversion in 
the steady-state. Furthermore, large refractive indices are feasible as well as the determination of 
the diagonal matrix elements via the absorption or emission spectra. The results apply to available 
biomolecular, spin or asymmetric quantum dot systems. 
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I. INTRODUCTION 

Population inversion in a few-energy-level quantum 
system is strongly linked with its optical response and 
has resulted in traditional masing or lasing effects be¬ 
ing successfully demonstrated PHI . Subsequently, enor¬ 
mous economical and technological progress was trig¬ 
gered due to quantum electronics. However at present, 
atomic steady-state population inversion is achievable es¬ 
sentially only in pumped two-level atomic systems involv¬ 
ing efforts via extra transitions with additional photon 
sources or environmental vacuum modifications |J|. Ad¬ 
ditionally, two-level quantum dot systems may also ex¬ 
hibit population inversion due to extra phonon induced 
decay rates Q. Meanwhile, lasing without population 
inversion operates as well with the help of induced quan¬ 
tum coherences f@|. Moreover, free-electron laser sources 
involving inversions of momentum states exist at higher 
frequencies leading to original effects |Tj. 

Over the past decade a range of experiments were per¬ 
formed demonstrating quantum effects in biological sam¬ 
ples [§}. Particularly, energy transfer through quantum 
coherence in photosynthetic systems [S] was observed in 
0, El, while long-living quantum coherences survive 
in biological complexes even under normal conditions at 
room temperature flU ]. Quantum coherence and entan¬ 
glement in the processes of mag neto-reception of the sur¬ 
rounding magnetic field fl3l.Il4l| were examined, too fl5j . 
The single molecule fluorescence spectroscopy and the 
emitted photon quantum statistics are further excellent 
tools for research in quantum biology m, S3- Further¬ 
more, coherent control of an effective two-level system in 
a non-Markovian biomolecular environment was investi¬ 
gated [h|. Additional quantum effects in biochemical 
systems are discussed, for instance, in [l9j. An impor- 
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tant issue raised in this context is to image tiny objects, 
such as biological cells or organic molecules and, there¬ 
fore, highly refractive biological media are required [20j. 

Inspired by these remarkable advances in quantum bio¬ 
chemistry, here, we put forward a novel scheme that en¬ 
ables the creation of population inversion in certain bio¬ 
logical samples acting as two-level systems with all possi¬ 
ble dipole matrix elements being nonzero. The effect oc¬ 
curs due to induced quantum coherences which arise from 
the difference of the permanent dipole moments in the ex¬ 
cited and ground states, respectively. This may allow for 
lasing or amplifying as well as optical switching devises in 
biomolecular materials. Large index of refractions with¬ 
out absorption are also feasible which may lead to an 
enhanced optical imaging resolution of the biomolecular 
sample due to a reduced probe-field wavelength inside the 
medium. Both the absorption or emission spectra can be 
used to extract the values of permanent dipoles. 

The article is organized as follows. In Sec. II we de¬ 
scribe the analytical approach and the system of interest, 
whereas in Sec. Ill we analyze the obtained results. A 
summary is given in Sec. IV. 

II. THEORETICAL FRAMEWORK 

We consider a two-level system possessing permanent 
dipoles and interacting with two external coherent laser 
fields. The first laser is near resonance with the tran¬ 
sition frequency of the two-level sample while the sec¬ 
ond one is close to resonance with the dressed-frequency 
splitting due to the first laser, respectively (see Figure 
[[]). As a concrete system, we may consider gamma¬ 
globulin macromolecules [19|, [2jJ with the transition fre¬ 
quency 0 J 21 ~ 4.8 x 10 15 Hz, transition dipole moment 
d ~ IDebye and the difference between the diagonal 
dipole moments in the upper and lower bare states given 
by 1^22 — dn| ~ lOODebye. However, the analytical for¬ 
malism applies equally to spin [22}, asymmetric semicon¬ 
ductor quantum dot p3j or other alternative [24| systems 
promising wider applications. The Hamiltonian describ- 
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FIG. 1: Energy diagram of a two-level emitter with nonzero 
values of all involved dipole matrix elements d al 3, {a, ft = 
1,2}. (a) A moderately intense laser field of frequency ujl 
resonantly interacts with the molecular sample leading to dy¬ 
namically Stark splitting of its energy levels, (b) A second 
coherent source of frequency co, close to generalized Rabi fre¬ 
quency due to first laser, is applied leading to transitions 
among the dressed-states. (c) The ’’double dressed-states” 
with 2 Gr being the corresponding Rabi splitting frequency. 


ing such a model, in a frame rotating at the first laser 
frequency wi, and in the dipole approximation is: 

H = Y, ~ u L )a\a k + HAS, + hfl(S+ + S~) 

k 

+ HGS Z cos ( cot ) + i ^2(gk ■ d){a\S~ - akS + ). (1) 

k 

In the Hamiltonian (JT]) the first three components are, re¬ 
spectively, the free energies of the environmental electro¬ 
magnetic vacuum modes and molecular subsystems to¬ 
gether with the laser-molecule interaction Hamiltonian. 
There, H = dEi/{2h) is the corresponding Rabi fre¬ 
quency with d = g ?2 1 = d ±-2 being the transition dipole 
moment while E\ is the amplitude of the first laser field. 
The fourth term accounts for the second laser interacting 
at frequency co and amplitude E 2 with the molecular sys¬ 
tem due to presence of permanent dipoles incorporated in 
G, i.e., G = )d 22 — dii)E 2 /h. The last term describes the 
interaction of the molecular subsystem with the environ¬ 
mental vacuum modes of the electromagnetic field reser¬ 
voir. Further, g = ^2nhj0k/Ve\ is the molecule-vacuum 
coupling strength with e\ being the photon polarization 
vector and A £ {1,2} whereas V is the quantization vol¬ 
ume. A = w 2 i — Wi is the laser field detuning from the 
molecular transition frequency co 2 \. The molecule bare- 
state operators 5 ,+ = |2)(11 and S~ = [<S ,+ ] 1 ' obey the 
commutation relations [S'" 1 ", S' - ] = 2 S z and [S^S 1 * 1 ] = 
±S ± . Here, S z = (|2)(2| — |l)(l|)/2 is the bare-state 
inversion operator. |2) and |1) are, respectively, the 
excited and ground state of the molecule while a/ and 
ak are the creation and the annihilation operator of the 
kth electromagnetic field mode, and satisfy the standard 
bosonic commutation relations, namely, [afc,aj_,] = 5 kk >, 
and [ak,a k '] = [a[.,aj],] = 0 [2514281 ]. Notice that the 


Hamiltonian m is incomplete. The following term 

H = hflS + e zulLt cos (cot) + HGS Z cos (wz,f)/2 + H.c., (2) 

will not be taken into account. Here, the first term de¬ 
scribes the interaction of the second laser with the molec¬ 
ular system due to the transition dipole d, while the sec¬ 
ond one represents the interaction of the first laser with 
the molecule because of permanent dipoles and, hence, 
H = dE 2 /h whereas G = ( d 22 — dn)Ei/H. Since we as¬ 
sume realistic conditions that tl <C wl±w [27|, [28| as well 
as {G, w} <C ul the Hamiltonian (O can be considered 
as fast oscillating and, therefore, neglectable. 

In the following, we shall consid er a regime w here the 
generalized Rabi frequency fl = ^/fl 2 + (A/2 ) 2 is larger 
than the single-molecule spontaneous decay rate as well 
as the coupling due to permanent dipoles, i.e. Q 7 and 
H > G. In this case it is more convenient to describe our 
system in the semi-classical laser-molecule dressed-state 
picture due to the first applied laser: 

|2) = cosd|2) — sindll) and |1) = cos dll) + sind|2), (3) 

with tan 29 = 2H/A. Here |2) and |1) are the correspond¬ 
ing upper and lower dressed states, respectively (see Fig¬ 
ure Eb). Applying the dressed-state transformation to 
the Hamiltonian (JTJ) one arrives at the following Hamil¬ 
tonian represented in a frame rotating also at the second 
laser field frequency, i.e. co, 

H = — 0JL)a{ak + hAR z — hG(R + + R ~) 

k 

+ i gk ■ d){a\{sva.29R z /2 + cos 2 9R~e~ lult 

k 

- sin 2 9R + e iuJt ) — iT.c.}, (4) 

where A = Cl — w(l — G 2 /co 2 )/2 = Cl — co/2 and G = 
(G/4)sin2d, and we have performed the rotating wave 
approximation with respect to co, i.e., we have assumed 
that u>G. Eliminating the vacuum modes of the elec¬ 
tromagnetic field reservoir in the usual wav by adopt¬ 
ing the Born-Markov approximations [251 - 1281 ] one arrives 
then at the following dressed-state master equation: 

p(t) + i[^ R z - G(R + + R~), p]=- : f sin 2 2 9[R z ,R z p\ 
— 7 + cos 4 9[R + , R~p] — 7 _ sin 4 9[R~, R + p] + H.c. (5) 

Here 70 ,± are the single-qubit spontaneous decay rates 
corresponding to dressed-state frequencies 00 l and uj, ± 
2f2, respectively. The new quasi-spin operators, i.e., 
R+ = 12}(T|, R~ = [i?+]t and R z = |2}(2| - |1)(1| are 
operating in the dressed-state picture. They obey the 
following commutation relations: [i? + ,i? - ] = R z and 
[R z , i? 41 ] = ±2i? ± . Notice that in the above master equa¬ 
tion, we have neglected the rapidly oscillating terms in 
the spontaneous emission part - an approximation valid 
when {H,w} > { 70 ,7±}- 

In the following section, we will discuss our results, i.e., 
the possibility to create bare-state population inversion 
as well as high refractive media. 
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FIG. 2: (color online) The steady-state dependence of the 
bare-state inversion operator {S z )/N versus the scaled param¬ 
eter A/(2fi). The solid line is plotted for G /7 = 0, the long- 
dashed one is for G /7 = 16 while the short-dashed curve cor¬ 
responds to G /7 = 24. Other parameters are: Q/(N , y) = 45 
and u/(N 7 ) = 100 with 7 ~ 2.6MHz as feasible for gamma¬ 
globulin. (a) IV = 1 whereas (b) TV = 50 while molecules are 
dense enough to allow for collectivity. 


III. RESULTS AND DISCUSSION 

One can observe from Eq. ([5]) that permanent dipoles 
lead to the appearance of a pumping term among the 
dressed states. This term contributes to a completely un¬ 
expected behavior of the two-level system. In particular, 
as we shall see, it can induce pumping the system into an 
inverted state. The system of equations for the dressed- 
state inversion and dressed-state polarization operators 
can be obtained from the master equation (|5|). namely, 

(Rz) = ~2iG((R-)-(R+))^2r + {R z ) + 2T_, 

< R + ) = (21A — T)(R + ) + iG(R z ), ( 6 ) 

with ( R~ ) = [(l? + )]b Here, T±= 7 (sin 4 0 ± cos 4 0), T = 
T+ +7 sin 2 28 and we have considered that 70 = 7 ± = 
7 which is the case for a free-space setup. The mean- 
value of the bare-state inversion operator (S z ) can be 
represented via dressed-state operators as follows: 

(S z ) = cos29(R z )/2 - sin2 9({R + ) + (R~))/2. (7) 

From the system of equations ([ 6 ]) one immediately ob¬ 
tains the steady-state relations: 

(R + ) = iG(R z )/(T - 2iA) and (R~) = [(i? + )] t . ( 8 ) 
Inserting ((5J) in Eq. m one arrives at: 

(S z ) = (cos 20 + 4AGsin20/[T 2 + (2A) 2 ])(7? z }/2, (9) 
where, again, from Eqs. m one has that: 

(R z ) = 2T_/[2T+ + (2G) 2 T/(T 2 + (2A) 2 )]. (10) 

An interesting result here is the non-zero value for the 
dressed-state coherences (R^, see Eq. ((SJ). In the ab¬ 
sence of permanent dipoles, i.e. G = 0, these quantities 
are zero in the moderately intense pumping regime con¬ 
sidered here. Therefore, new physics is expected due to 
existence of permanent dipoles in bichromatic pumping 


fields. In particular, Figure [2j a) shows the mean-value of 
the single-molecule bare-state inversion operator for par¬ 
ticular parameters of interest. Steady-state inversion in 
the bare states is achieved when (S z ) > 0 and it occurs in 
the presence of permanent dipoles. On the other hand, 
Figure[3Ka) depicts the real part of the mean-value of the 
dressed-state coherence operator ( R + } in steady-state. 
The minima observed in these behaviors correspond to 
an inverted molecular bare-state system (compare Fig.[2ji 
and Fig.[3Jr). Thus, inversion occurs due to the real part 
of the dressed-state coherences which can be nonzero in 
our system. 

One can apply the double dressed-state formalism (see 
Figure |T|) in order to obtain further information on our 
system, namely, 

12 ) = cos</>| 2 ) + sin 0 | 1 ), | 1 ) = cos</>|l) — sin 0 | 2 ). ( 11 ) 

This approach is particularly useful to account for vac¬ 
uum induced collective effects among the emitters and 
their corresponding influences on molecular dynamics. 
Introducing (fTTl) in the Hamiltonian (|1| and, again, elim¬ 
inating the degrees of freedom related with the environ¬ 
mental vacuum modes in the Born-Markov approxima¬ 
tions one arrives at the double dressed master equation: 

p{t) + iG R [R z ,p] =-T 0 {R z ,R z p]-T + [R + ,R~p\ 

- t-[R~,R + p]+ H.c. (12) 

Here, f 0 = 7 ( 07 ,) sin 2 20 cos 2 20/4 + sin 2 20 ( 7(07 + 
w)cos 4 0 + 7 (wl — w)sin 4 0}/4, f + = 7 (07 + 

2 Gr) sin 2 2(f) sin 2 20/4 + 7 ( 1 x 7 + oj + 2 Gr) cos 4 0 cos 4 0 + 
7(0 7 — o> + 2 G r ) sin 4 0 sin 4 0 and T_ = 7(07 — 

2G r ) sin 2 2/> sin 2 20/4 + 7(07 + w — 2G',r) cos 4 0 sin 4 0 + 
7 (07 — o> — 2 G r ) sin 4 0 cos 4 0 with cot 20 = A/G, and 
Gr = VA 2 + G 2 . The new operators, i.e., R + = 
|2)(1|, R~ = [I? + ]t and R z = |2)(2| — |1)(1| are op¬ 
erating in the double dressed-state picture obeying the 
following commutation relations: = R z and 

[Rz.lG 1 ] = ±2I? ± . The master equation (fT2l) contains 
only slowly varying terms in the spontaneous emission 
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FIG. 3: (color online) The single-molecule steady-state behav¬ 
iors of (a) the real part of the dressed-state coherences ( R + ) 
and (b) the double dressed-state inversion operator ( R z ) ver¬ 
sus the scaled detuning A/(212). The solid line is plotted for 
G /7 = 16 while the dashed curve corresponds to G /7 = 24. 
Here fi /7 = 45 and 0/7 = 100. 
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damping, that is, we have assumed that Gr 7 (w), with 
7 (A) = 2d 2 £j 3 /(3hc 3 ) being the single-molecule sponta¬ 
neous decay rate corresponding to the double dressed- 
state frequency w (see Figure [1]). 

The steady-state solution of Eq. (fl 2 l) can be chosen in 
the form: 


p = Z 1 exp[-r]R z ], (13) 

where the normalization Z is determined from the re¬ 
lation Tr{p} = 1. Inserting (1131) into (1121) one obtains 
77 = ln(f+/r_)/ 2 . Using the relations © and (fill) we 
arrive at the following expression for the mean value of 
the bare-state inversion operator (S z ) represented via 
the double dressed-state inversion operator (R z ), respec¬ 
tively, 

(S z )=cos[2(d-cj ) )}(R z )/2. (14) 


The steady-state expression for the double dressed-state 
inversion operator can be obtained with the help of 
Eq. (HI and the coherent molecular state |s) = \N — s, s) 
which denotes a symmetrized N —molecule state in which 
N — s molecules are in the lower double dressed-state |1) 
and s molecules are excited to the upper double dressed- 
state 12), respectively [23] ■ Thus, 


(Rz) 


A (x 1+N — !) — (! + N)(x — 1) 
2 - 1 (s-l)(a: 1 + JV -l) 


(15) 


where x = T + /T_. In this case the molecular system 
has to occupy a volume with linear dimensions of the or¬ 
der of the smallest relevant emission wavelength or less. 
However, this restriction may be relaxed for certain ge¬ 
ometries The molecular operators entering in 

Eqs. (I12II15D are then collective ones, i.e. R-f 

as well as R z = R*j- Figure [21b) shows the 

bare-state inversion operator ( S z )/N , based on the ex¬ 
pressions m and m for a collection of N = 50 
molecules in a volume (2 ttc/u) 2 i ) 3 or a molecular den¬ 
sity N ss 7.8 x 10 14 cm~ 3 . Almost full inversion in 
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FIG. 4: (color online) The steady-state dependence of the 
linear susceptibility x( v p) [in units of Nd 2 /(Try)] versus scaled 
detuning Ap/ 7 . The solid black curve corresponds to the 
imaginary part (absorption spectrum) while the long-dashed 
blue line to the real part of the susceptibility, respectively, 
(a) A/(2fi) = 0.43 whereas (b) A/(2fi) = —0.43. Other 
parameters are: fi /7 = 45, w /7 = 100, Gj 7 = 16. 


the steady-state is achieved when we consider that the 
double-dressed decay rates 7 (w) are equal and denoted 
by 7 . Furthermore, abrupt population behaviors are ob¬ 
served as well which may be used for engineering opti¬ 
cal switching devices with switching times of the order 
of (Aby) -1 . Note that the mean values of non-diagonal 
terms resulting from Eq. m are zero in the steady-state. 

We proceed by calculating the refractive properties of 
a very weak held probing the strongly driven molecular 
sample around the first laser’s frequency u>l- The lin¬ 
ear susceptibility \(v) of the probe held, at frequency 

v, can be represented in terms of the Fourier transform 
of the average value of the two-time commutator of the 
molecular operator as 

pOO 

X (v) = j:Nd 2 J o dre^-^([S-(T),S+]) s . (16) 

Note that the subindex s means steady-state. Inserting 
the dressed-state transformations © and m in m, 
in a frame rotating also at the second laser’s frequency 

w, and using the master equation (1121) one arrives at the 
following expression for the susceptibility, namely, 

xM = ^Nd 2 (R z ){^sm 2 26sm 2 2</>(xi(A p , 2Gr) 

- Xi(A p , - 2 G r )) + cos 4 d(sin 4 <fo 2 (A p , 2 G R , -w) 

— cos 4 (/>X 2 (A p , — 2 Gr, —ijo)) + sin 4 0 (cos 4 <j) 

x X 3 (Ap, 2Gr,lu) — sin 4 1 ^x 3 (A p , —2Gr, w))}, (17) 

where xi(A p ,a:) = (f s + *(A P + :r))/Jf 2 + (A p + x) 2 ], 
X 2 (A p ,x,-y) = (r s +i(Ap + x-y))/[r2 + (Ap + x-y) 2 ] 
and x 3 (Ap, x, y) = (r s +*(Ap+j:+y))/[r 2 -j-(Ap+x-t-?/) 2 ] 
whereas A p = v — ujl while T s = 4Fo + r + + T_. Fig¬ 
ure a shows the steady-state behavior of the linear sus¬ 
ceptibility when the molecular sample is probed with a 
weak coherent held of frequency v. Both, positive or 
negative dispersions without absorption are clearly vis¬ 
ible around Ap /7 = 0 which may lead to enhanced or 
reduced refractive indices applicable for optical imaging, 
lithography or negative refraction processes in dense me¬ 
dia [29j . In particular, the index of refraction close to 
vanishing absorption n{y) ~ y/l + x takes values n > 2 
for N = 10 17 cm -3 and A/(2$7) = —0.43. These depen¬ 
dences for the susceptibility xi v ) are easily understood 
in the double dressed-state picture (see Fig. m ■ Partic¬ 
ularly, for A/(2f2) = 0.43 there are more molecules in 
the upper double dressed state |2), i.e. (R z ) > 0 (see 
Fig. ©3). As a consequence, at frequency wl_+w + 2G,r 
one has gain while at frequency ujl — to — 2 Gr we have 
absorption, see the lateral dip/peak in Figure © 1 ). Sim¬ 
ilarly one can explain the whole structure shown in Fig¬ 
ure ©a, b). An interesting issue here is the structure at 

± 2 Gr which may help to extract the value of perma¬ 
nent dipoles, i.e. |o ?22 — dn|. If one inspects the absorp¬ 
tion spectrum shown in Fig. ©then the frequency sepa¬ 
ration between the first maximum and minimum around 
Ap = 0 is equal to 4GA. Gr involves the difference of 
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permanent dipoles |c ?22 — rfn | • This, in principle, allows to 
determine |d 2 2 — dn| if all other involved parameters are 
known. Finally, the elastic photon scattering spectrum 
consists of three lines at {wl, ojl ± w}. The inelastic one 
may contain up to nine spectral lines in strict concor¬ 
dance with the double dressed-state formalism schemat¬ 
ically shown in Fig. ©)• Suppression of a spectral line 
at the frequency of the strongly driving laser also oc¬ 
curs f30j . As well, the emission spectrum permits de¬ 
termination of the diagonal dipole matrix elements (see 
also Ref. 0). The magnitudes of permanent dipoles are 
required for e.g. interpretation of biological images as 
well as for biological high-harmonic or ultrashort pulse 
generation processes via laser pumped media possessing 
permanent dipoles (32 (see also [33]). 

IV. SUMMARY 

Summarizing, we have investigated the steady-state 
quantum dynamics of laser pumped two-level molecu¬ 
lar samples with broken inversion symmetry. We have 
demonstrated population inversion in the bare states due 


to induced coherences which in turn depend on the mag¬ 
nitude of permanent dipoles. Vacuum mediated collec¬ 
tive effects among the two-level emitters considerably en¬ 
hance the molecular inversion. The values of the perma¬ 
nent dipoles can be inferred from the emission or ab¬ 
sorption spectra. Furthermore, the investigated system 
exhibits large positive or negative dispersion without ab¬ 
sorption facilitating applications including optical imag¬ 
ing, lithography and negative indices of refraction. The 
results apply especially to biomolecular, spin or asym¬ 
metrical quantum dot systems. 
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